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Abstract 



qq , Hypersensitivity to perturbation is a criterion for chaos based on the ques- 

tion of how much information about a perturbing environment is needed to 
keep the entropy of a Hamiltonian system from increasing. We demonstrate 
numerically that hypersensitivity to perturbation is present in the following 
quantum maps: the quantum kicked top, the quantum baker's map, the quan- 
. turn lazy baker's map, and the quantum sawtooth and cat maps. 
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I. INTRODUCTION 



Present research on quantum chaos focuses mainly on what Berry |IJ has called quantum 
chaology, namely "the study of semiclassical, but nonclassical, behavior characteristic of 
systems whose classical motion exhibits chaos." We pursue a different approach, which 
seeks to identify intrinsically chaotic features within quantum dynamics. A straightforward 
application of this approach, which attempts to generalize the classical notion of sensitivity 
to initial conditions, fails immediately, for two reasons: first, there is no quantum analogue 
of the classical phase-space trajectories in terms of which sensitivity to initial conditions 
is defined and, second, the unitarity of linear Schrodinger evolution precludes sensitivity 
to initial conditions in the quantum dynamics of state vectors. A classical Hamiltonian 
system has, however, in addition to its description in terms of trajectories governed by the 
Hamilton equations of motion, an equivalent description in terms of Liouville probability 
densities governed by the Liouville equation. Because it conserves phase-space volumes, 
the linear Liouville equation shares with the Schrodinger equation a lack of sensitivity to 
initial conditions. Again, in Berry's words ||, "the fact that this equation is linear and 
yet preserves the chaos of the trajectories disposes of the argument that the linearity of 
the Schrodinger equation is responsible for the absence of chaos in quantum mechanics." 
We have recently established a criterion for classical chaos in terms of Liouville probability 
densities ||-|5[] , which we call hypersensitivity to perturbation and which has a straightforward 
quantum analogue @,§,[7|- 

Hypersensitivity to perturbation for a quantum system is defined in purely quantum 
terms, not in the semiclassical domain. Though quantum systems show no sensitivity to 
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initial conditions, due to unitarity, they do show what one might call sensitivity to param- 
eters in the Hamiltonian, as has been demonstrated for the quantum kicked top by Peres 
H|9|]. Our approach to quantum chaos could be viewed as a generalization of Peres's work: 
whereas Peres compares the time evolution of the same initial Hilbert-space vector for two 
slightly different parameter values in the kicked-top Hamiltonian, we study the distribution 
of Hilbert-space vectors arising from the many possible ways in which the Hamiltonian time 
evolution can be perturbed through interaction with an environment. Our approach, how- 
ever, is embedded in a broader information-theoretic framework 0,|5],|7| that is quite different 
from the viewpoint expressed in ||. 

In this paper we present numerical evidence that hypersensitivity to perturbation is 
present in a variety of quantum maps. We make no attempt to develop a general theory of 
quantum hypersensitivity to perturbation; rather, our attitude in this paper is exploratory: 
by investigating numerically the perturbed dynamics of several quantum systems, we seek 
to determine whether hypersensitivity to perturbation provides a criterion for distinguish- 
ing regular from chaotic dynamics in quantum systems. The paper is organized as follows. 
Section [TT] defines hypersensitivity to perturbation and reviews the main results for clas- 
sical Hamiltonian systems. Section |T| introduces the model of the quantum kicked top 
interacting with a perturbing environment. Section |1V A| explains how hypersensitivity to 
perturbation can be related to the distribution of vectors in Hilbert space; Sec. [IV B| contains 
numerical results for the quantum kicked top. In Sec. [V| we present numerical evidence for 
hypersensitivity to perturbation in a variety of quantum maps on the unit square. Section [V] 
concludes with brief discussion. 



II. HYPERSENSITIVITY TO PERTURBATION 

Hypersensitivity to perturbation, in either classical or quantum mechanics, is defined 
in terms of information and entropy. The entropy H of an isolated physical system (Gibbs 
entropy for a classical system, von Neumann entropy for a quantum system) does not change 
under Hamiltonian time evolution. If the time evolution of the system is perturbed through 
interaction with an incompletely known environment S, however, averaging over the per- 
turbation typically leads to an entropy increase AHs- The increase of the system entropy 
can be limited to an amount AH to \, the tolerable entropy increase, by obtaining, from the 
environment, information about the perturbation. We denote by A/ min the minimum infor- 
mation about the environment needed, on the average, to keep the system entropy below 
the tolerable level AH to \. A formal definition of the quantities AH$, AH to \, and AJ min is 
given below. 

Entropy and information acquire physical content in the presence of a heat reservoir at 
temperature T. If all energy in the form of heat is ultimately exchanged with the heat 
reservoir, then each bit of entropy, i.e., each bit of missing information about the system 
state, reduces by the amount ksT In 2 the energy that can be extracted from the system in 
the form of useful work. The connection between acquired information and work is provided 
by Landauer's principle p^TTl , according to which not only each bit of missing information, 
but also each bit of acquired information, has a free-energy cost of /csTln2. This cost, 
the Landauer erasure cost, is payed when the acquired information is erased. Acquired 
information can be quantified by algorithmic information |T2|-|T^]. 
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We are now in a position to define hypersensitivity to perturbation. We say that a system 
is hypersensitive to perturbation if the information A/ min required to reduce the system 
entropy from AH S to AH to \ is large compared to the entropy reduction AH S — AH to i, i.e., 

A/ " in »1. (2.1) 



AHq - AH 



tol 



The information A/ min purchases a reduction AH$ — AH to \ in system entropy, which is 
equivalent to an increase in the useful work that can be extracted from the system; hyper- 
sensitivity to perturbation means that the Landauer erasure cost of the information is much 
larger than the increase in available work. 

Hypersensitivity to perturbation means that the inequality (|2.1| ) holds for almost all 
values of AH to \. The inequality ( |2.1| ) tends always to hold, however, for sufficiently small 
values of AH to \. The reason is that for these small values of AH to \, one is gathering enough 
information from the perturbing environment to track a particular system state whose en- 
tropy is nearly equal to the initial system entropy. In other words, one is essentially tracking 
a particular realization of the perturbation among all possible realizations. Thus, for small 
values of AH to \, the information A/ min becomes a property of the perturbation; it is the 
information needed to specify a particular realization of the perturbation. The important 
regime for assessing hypersensitivity to perturbation is where AH to \ is fairly close to AHs, 
and it is in this regime that one can hope that AJ min reveals something about the system 
dynamics, rather than properties of the perturbation. 

In || we have shown that a large class of classically chaotic systems shows exponential 
hypersensitivity to perturbation, which is characterized by an exponential growth with time 
of the ratio between information and entropy reduction, for a broad range of values of 
AH to \ below AHs- The exponential rate of growth is determined by the Kolmogorov- Sinai 
entropy of the classical dynamics. In preliminary work [[^,[7] we have applied the notion of 
hypersensitivity to perturbation to quantum systems and, through numerical simulations, 
have shown that hypersensitivity to perturbation seems to distinguish regular from chaotic 
quantum dynamics. The present paper contains a more detailed treatment of hypersen- 
sitivity to perturbation, including numerical evidence for a range of prototypic quantum 
systems. 

To define the quantities AHs, AH to \, and A/ min rigorously in the quantum case, consider 
a quantum system S coupled to an environment S. We assume that the total system S <S> £ 
is initially in a product state 

Ptotal(* = 0) = \^ )(M®P£ , (2.2) 

combining a pure system state \ipo) with a mixed environment state ps- The initial von 
Neumann entropy H of the system is zero. After evolving for a time t under the action of 
the total Hamiltonian H totah the total system is described by a density operator ptotai(^); 
which is generally not a product state. The state ps of the system S alone at time t is then 
obtained by tracing over the environment, 

As = tr £ (p t otai(t)) • (2.3) 
The increase in system entropy due to tracing over the environment is given by 
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A# 5 = -tr 5 (p 5 log 2 p s ) . (2.4) 

Now suppose a measurement is performed on the environment £ . The most general 
measurement is described by a positive-operator- valued measure (POVM) @|I7|]. We assume 
for convenience that the measurement is described by a POVM that has a discrete set of 
possible outcomes r with probabilities p r . For each measurement outcome r, we denote by p r 
the system state conditional on the measurement outcome r; the conditional system states 
obey the relation (for a justification, see Appendix A) 

^2Prf>r=PS- (2.5) 

r 

We further define the system entropy conditional on measurement outcome r, 

AH r = -tr s (p r log 2 p r ) , (2.6) 
the average conditional entropy 

AH = Y,PrAH r , (2.7) 

r 

and the average information 

AI = -$> r log 2 p r . (2.8) 

r 

Actually AI is a lower bound to the average algorithmic information needed to specify 
the measurement outcome r. One of us has shown |TJ|, however, that the minimum aver- 
age algorithmic information can be bounded above by AI + 1, so that AI is an excellent 
approximation to the minimum average algorithmic information. 

The minimum amount of information about the perturbing environment to keep the 
system entropy from increasing by more than the tolerable amount AH to \ can now be defined 
as 

where the infimum is taken over all environment measurements for which the average condi- 
tional entropy increase AH does not exceed AH to \. In other words, A/ min is the information 
about the environment it takes to reduce the entropy increase of the system from AH$ (the 
increase due to averaging over the perturbation) down to AH to i', i.e., AJ min is the minimum 
information about the environment needed to purchase a reduction AH$ — AH to \ in system 
entropy. 



III. QUANTUM KICKED TOP 
A. Kicked-top Hamiltonian 



The quantum model of the kicked top | JT8"| , |iT| describes a spin- J particle — i.e., an angular 
momentum vector hJ = h(J x , J y , J z ), where [ Jj, Jj] = itijkJk — whose dynamics in (2J + 1)- 
dimensional Hilbert space is governed by the Hamiltonian 
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+oo 

H top (t) = (hp/T)J z + (hk/2J)J 2 x £ 5(t-nT) 

n=— oo 



(3.1) 



The free precession of the spin around the z axis (first term in the Hamiltonian) is interrupted 
periodically by sudden kicks or twists at times nT with twist parameter k (second term in 
the Hamiltonian). The angle of free precession between kicks is given by p. In this paper 
we use p = 7r/2. 

We look at the time evolution of an initial Hilbert-space vector |^ ) at discrete times 
nT. After n time steps, the evolved vector is given by 

|Vn> = f; op |^ ) , (3-2) 

where T top is the unitary Floquet operator 

f top = e - l ( fc / 2J ) J 'e-^/ 2 . (3.3) 

The model has a conserved parity, which for half-integer J takes the form 

S = -ie~ inX (3.4) 

and which permits factorization of the matrix representation of the operator T top into two 
blocks. Starting from a state with definite parity, the whole dynamical evolution occurs 
in the invariant Hilbert subspace with the given parity. For half-integer J, the dimension 
of the even-parity subspace (eigenspace of S with eigenvalue 1) is J + |. In this paper, 
we work with J = 511.5 in the 512-dimensional even-parity subspace or with J = 63.5 
in the 64-dimensional even-parity subspace. We consider only the projection of the initial 
vector in the even subspace. Numerical evidence and symmetry considerations || suggest 
that no additional insight is gained by including the odd-parity subspace. In any case, 
the restricted model can be regarded as a quantum map in its own right, which can be 
investigated independently of the behavior of the complete kicked-top model. 



B. Choice of initial vector 



Depending on the initial condition, the classical map corresponding to the Floquet op- 
erator ( |3.3| ) displays regular as well as chaotic behavior [|19|] . Following M , we choose initial 
Hilbert-space vectors for the quantum evolution that correspond to classical initial condi- 
tions located in regular and chaotic regions of the classical dynamics, respectively. 

The classical map corresponding to the quantum map (|3.3|) is obtained by introducing 
the unit vector (X, Y, Z) = J/J and performing the limit J — > oo. One obtains JT9| 



X' = -Y , 

Y' = XcoskY + ZsinkY , (3.5) 
Z' = Z cos kY - X sin kY , 

which is an area-preserving map of the unit sphere, i.e., an area-preserving map on the 
configuration space of a classical spin with fixed magnitude. Depending on the value of the 
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twist parameter k, this map has regions of chaotic behavior interspersed with regular regions 
associated with elliptic cyclic points. In this paper, we are interested in two cyclic points of 
the map for k = 3. One is an elliptic fixed point of period 1, located at 



Z = cos 9 = 0.455719 , p = 3vr/4 , (3.6) 

where we use spherical coordinates 9 = cos _1 Z and tp = ta.n~ 1 (Y/X). For k = 3 the 
elliptic fixed point ( fT6| ) is surrounded by an oval-shaped regular region of regular dynamics, 
extending about 0.3 radians in the (^-direction and about 0.5 radians in the ^-direction. 
The other cyclic point of interest to us is a hyperbolic fixed point of period 4, which has a 
positive Lyapunov exponent. It is located at fl20| 



Z = cos 9 = 0.615950 , <p = vr/4 , (3.7) 

in the middle of a chaotic region 

We want to choose initial vectors for the quantum evolution that correspond as closely as 



possible to the classical directions (|3.6|) and ( |3.7| ). For this purpose, we use coherent states 



\9,<p) [PTJ-|23|, defined by the relation 

n-J\6,(p) = J\6,(p), (3.8) 

where n is the unit vector pointing in the direction given by 9 and p. All coherent states 
can be generated by an appropriate rotation of the state \J,M — J) = \9 = 0, p = 0), 
where \J, M) (M = — J, . . . , J) is the simultaneous eigenstate of J 2 and J z with eigenvalues 
J( J + 1) and M, respectively. 

The distance between two normalized vectors \ipi) and \ip 2 ) is defined as the Hilbert-space 
angle 

S (|^ 1 ),|^ 2 )) = cos- 1 (|^i|^ 2 )|) (3.9) 



between the two vectors |24]. Consider two coherent states \9,<p) and \9',p'}. In terms of 
the angle a between the directions (9, p) and (9', <p'), the distance between the two coherent 
states is given by |25| 



cos[s(|0,p>,|0V»] = \(9,p\9',p')\ = [cos(a/2)] 2J ~ exp(- Ja 2 /A) , (3.10) 

where the approximation is valid for large J. Two coherent states can therefore be regarded 
as roughly orthogonal if a > [H]. The size of the coherent state \9, p) is conveniently 

defined in terms of the Q function 

Q e ^(9',p') = \(9',p'\9,p)\ 2 = [cos(«/2)] 4J = Q(a) . (3.11) 

Since Q(2J~ 1 ^ 2 ) ~ e~ 2 Q(0), the Q function of the coherent state \9, p) is very small outside 
a region of radius centered at the direction (9,<p). For the value J = 511.5 used in 

this paper to study regular behavior, one finds a radius of 2J -1 / 2 ~ 0.09 radians, less than 
the size of the regular region around the elliptic fixed point 
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C. Interaction with an environment 



Suppose now that the kicked top interacts with an environment consisting of a collection 
of degenerate two-state systems. For the interaction Hamiltonian we assume the form 

+00 

#int = hgj z fo)n6(t -nT-e). (3.12) 

n=— 00 

Here (a^) n is a Pauli operator for the nth two-state system, g is the interaction strength, and 
e is an infinitesimal time interval. The kicked top thus interacts with a single environment 
system at a time, interacting with the nth two-state system just after the nth twist, i.e., 
just after the nth unperturbed time step. The nth perturbed time step includes the nth 
unperturbed time step and the immediately subsequent coupling to the nth two-state system. 
The total Hamiltonian for system and environment is 

-f^total — H top + if int , (3.13) 

where H top is given by Eq. ( |3.1[) with p = tt/2. 

As was explained in Sec. [TJ, we take the initial state of the total system to be of the form 
|'0o}('0o| ® Pe, where here the initial state \ipo) of the system is a coherent state centered at 
one of the fixed points (|3.6|) and ( |3.7|) and the initial state of the environment, for the n 
two-state systems that are relevant for the first n time steps, is 

h= E 2- n \h)---\l n )(l n \---(h\ . (3.14) 
Ji ,...,*„£ {-l.+l} 

Each of the environment oscillators is with equal probability in the "lower" state | — 1), where 
a 3 \— 1) = — |— 1), or in the "upper" state |+1), where cr 3 |+l) = | + 1). 

Now suppose the kicked top and the environment evolve together for n time steps. We 
want to make a measurement on the environment in order to reduce the entropy increase of 
the system. The form of the interaction Hamiltonian ( |3.12|) suggests that the natural basis 



for a measurement of the environment state is the simultaneous eigenbasis of the operators 
(a^) m . If the environment is measured in this basis and is found to be in the state • ■ ■ \l n ) 
(l m G {— 1, +1}, m = l,...,n), then the system state after n perturbed steps is given by 

Wn) = T top (l n )---f top (h)\ip ) , (3.15) 

where 

f top {lm) = e - i9lmJ >e- l(k/2J)J *e- lwX/2 = e~ i9lmX f top (3.16) 

is the unperturbed Floquet operator ( |3.3| ) followed by an additional rotation about the z 
axis by angle gl m . Notice that T top (/ m ) commutes with parity ( |3.4| ) and, hence, does not 
couple odd- and even-parity subspaces. 

The effect of the coupling to the environment, together with the measurement, is to 
produce a binary stochastic perturbation at the end of each time step. The top rotates by 
+g if the two-state system is found in the upper state and by —g if it is found in the lower 
state, the two possibilities occurring with equal probability. After n steps, there are 2 n 
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different equally probable measurement results, each leading to a pure system state. Since 
there is no simple commutation relation between the operators T top (— 1) and T top (+l), these 
2 n system states are all different. Any measurement that fails to distinguish completely 
between different environment states • • • |/ n ) leads therefore to a mixed system state. To 
keep the system state pure, one has to specify which of the 2 n equally probable alternatives 
is realized. The minimal information about the environment needed to keep the system 
entropy at the level Aif t oi = is therefore A/ m i n = n bits. This means that the information 
about the environment needed to keep track of the exact system state is unbounded as the 
number of steps increases, an effect that is independent of the exact nature of the system 
dynamics. 

The value of A/ min for AH to \ = 0, therefore, is not a useful measure of the difficulty in 
keeping the system entropy small. We must determine how fast A/ min decreases if AH to \ is 
increased. Unfortunately, it is difficult to find an optimal measurement leading to minimal 
information for an arbitrary value of AH to \. It is possible, however, to find a nearly optimal 
measurement based on the following idea. The interaction with the environment Q3.12Q is 
of a special form that can be described in terms of a randomly perturbed unitary evolution 
operator. Each of the 2 n possible measurement results corresponds to one of 2 n different 
perturbation histories obtained by applying every possible sequence of perturbed unitary 
evolution operators T top (— 1) and T top (+l) for n steps. Consider the list C of 2 n vectors 
generated by all perturbation histories. In the next section, we show how a nearly optimal 
measurement can be found in the class of measurements whose outcomes partition the list 
C into groups of vectors. 

IV. DISTRIBUTION OF VECTORS IN HILBERT SPACE 

From this section on, we always assume that the interaction with the environment can 
be described by a stochastic unitary time evolution. In the first part of this section, we 
explore the relation between hypersensitivity to perturbation and the distribution of the 
vectors in Hilbert space resulting from different realizations of the stochastic evolution. In 
the second part, we use this relation to investigate hypersensitivity to perturbation in the 
quantum kicked top. 

A. General theory 

Given the assumption of a stochastic unitary time evolution, the interaction with the 
environment at a given time can be described by a list £ = ■ ■ ■ , \iPn)) of N vectors in 

ZJ-dimensional Hilbert space, together with probabilities qi, ■ ■ ■ ,qN, each vector in the list 
corresponding to a particular realization of the perturbation, which we call a perturbation 
history. Averaging over the perturbing environment leads to a system density operator 

N 

As = E^#j>(^l , (4-1) 

with entropy AH$ as given by Eq. ( |2.4|) . Consider the class of measurements on the envi- 
ronment whose outcomes partition the list L into R groups, labeled by r = 1, . . . , R. We 
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denote by N r the number of vectors in the rth group 

(Er=i N r = N). The N r vectors in the 
rth group and their probabilities are denoted by \ip[), ■ ■ ■ , \ij) r N ) and q[, . . . , q r N , respectively 
The measurement outcome r, occurring with probability 

N r 

Pr = Etf. ( 4 - 2 ) 

i=l 

indicates that the system state is in the rth group. The system state conditional on the 
measurement outcome r is described by the density operator 

N r 

pr^p^E^mm- (4.3) 

i=l 



Using Eqs. ( f4.2|) and (^4.3|) , we define the conditional system entropy AH r , the average 
conditional entropy AH, and the average information AI as in Eqs. ( |2.6| )-( |2~"5| ). 

We now describe nearly optimal measurements, i.e., nearly optimal groupings, in the 
sense of the discussion in Sec. |IT| Given a tolerable entropy AH to \, we want to partition 
the list of vectors C into groups so as to minimize the information AI without violating 
the condition AH < AH to \. To minimize AI, it is favorable to make the groups as large 
as possible. Furthermore, to reduce the contribution to AH of a group containing a given 
number of vectors, it is favorable to choose vectors that are as close together as possible 
in Hilbert space. Here the distance between two vectors and \ip2) can be quantified in 
terms of the Hilbert-space angle (|3.9|). Consequently, to find a nearly optimal grouping, we 
choose an arbitrary resolution angle (0 < < vr/2) and group together vectors that are less 
than an angle apart. More precisely, groups are formed in the following way. Starting with 
the first vector |-0i) in the list C, the first group is formed of and all vectors in £ that 
are within an angle of The same procedure is repeated with the remaining vectors 
to form the second group, then the third group, continuing until no ungrouped vectors are 
left. This grouping of vectors corresponds to a partial averaging over the perturbations. 
To describe a vector at resolution level amounts to averaging over those details of the 
perturbation that do not change the final vector by more than an angle 0. 

For each resolution angle 0, the grouping procedure described above defines an average 
conditional entropy AH = Aif(0) and an average information AI = A/(0). If we choose, 
for a given 0, the tolerable entropy AH to \ = AH(<p), then to a good approximation, the 
information AJ m ; n is given by A/ m i n ~ A/(0). By determining the entropy Aif (0) and the 
information A/(0) as functions of the resolution angle 0, there emerges a rather detailed 
picture of how the vectors are distributed in Hilbert space. If A/(0) is plotted as a function 
of Aif (0) by eliminating the angle 0, one obtains a good approximation to the functional 
relationship between A/ min and AH to \. 

It is easy to see that information and entropy obey the inequalities 

A/(0) > A/(0) > AI(n/2) = , (4.4) 
= AH(0) < AH((f)) < AH(n/2) . (4.5) 



The first of the information inequalities in Eq. ( |4.4| ) follows from the fact that any group at 
resolution is the union of groups at resolution = 0; in words, the average information 
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needed to specify a group at resolution = is equal to the average information needed to 
specify a group at resolution plus the average information needed to specify = groups 
within the groups at resolution 0. The latter of the entropy inequalities in Eq. ( |4.5| ) is a 
consequence of the concavity of the von Neumann entropy f26|.|27j . A general theorem about 
average density operators p6||2~T|| shows that, for all 0, 



A/(0) + AH((p) > AH{n/2) . (4.6) 

In general, AJ(0) is a decreasing function of 0, whereas Ai?(0) is increasing. This mono- 
tonicity can sometimes be violated, however, because of discontinuous changes in the group- 
ing of vectors. 

As a further characterization of our list of vectors, we calculate the distribution g(<ft) 
of Hilbert-space angles = s(\ip), = cos _1 (| (ip\ip') |) between all pairs of vectors 
and l^'). For vectors distributed randomly in D-dimensional Hilbert space, the distribution 
function g(0) is given by [|7| 

9(4) = dVD [f )/d<P = 2(D - l)(sin0)— cos0 . (4.7) 

Here V^(0) = (sin0) 2 ^ D_1 ^V£) is the volume contained within a sphere of radius in D- 
dimensional Hilbert space, and Vd = nD ^ 1 /(D — 1)! is the total volume of Hilbert space. 
The maximum of g(4>) is located at = cos -1 (l/^/2p - 1) ); for large-dimensional Hilbert 
spaces, g(0) is very strongly peaked near the maximum, which is located at ~ n/2 — 
1 /y/2D, very near tt/2. 

As yet a further characterization of the way the vectors are distributed in Hilbert space, 
we want to define a quantity that indicates how many dimensions of Hilbert space are ex- 
plored by the vectors in the list C. A possible measure of the number of explored dimensions 
is the exponential of the entropy, 2 AHs . This quantity is bounded above by D, the dimension 
of Hilbert space, and gets smaller if the dimensions are occupied with different weights. For 
example, if two eigenvalues of p$ are close to 1/2 and all the others are close to zero, then 
2 AHs ~ 2, indicating that the vectors are essentially confined to a two-dimensional subspace. 
Unfortunately, if all the vectors in £ are confined to a small sphere with radius 7r/2, 
AH$ is necessarily small just because all the vectors in the group lie along roughly the same 
direction in Hilbert space; this is true even if the orthogonal components of the vectors 
are evenly distributed over all the orthogonal directions in Hilbert space. For example, 
the density operator describing a uniform distribution of vectors within a sphere of radius 
<C 7i / 2 has one dominating eigenvalue close to 1 and D — 1 eigenvalues that are all equal 
and close to zero 0. Clearly, in this case 2 AHs is not an adequate measure of the number 
of dimensions explored. On the other hand, if one disregarded the largest eigenvalue in this 
example, then the exponential of the entropy would become a useful measure of the number 
of explored dimensions. We therefore introduce the spread AH 2 as the entropy calculated 
with the largest eigenvalue of ps omitted. The spread is defined as 

where Ai > A 2 > . . . > are the eigenvalues of the density operator p s . By giving different 
weight to dimensions corresponding to different eigenvalues of ps, the quantity [2 Ai?2 ] turns 
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out to be a good indicator of the number of Hilbert-space dimensions explored by the vectors 
in the list C, independent of the size of the region occupied by the vectors ( \x] denotes the 
smallest integer greater than or equal to x). In our analysis of the numerical results, we 
identify the number of dimensions explored by the entire set of N vectors with the integer 
n d =\2 AH ^. 

To investigate if a quantum map shows hypersensitivity to perturbation, we use the 
following numerical method. We first compute a list of vectors corresponding to different 
perturbation histories. Then, for about 50 values of the angle </> ranging from to tt/2, we 
group the vectors in the nearly optimal way described above. Finally, for each grouping and 
thus for each chosen angle 0, we compute the information AI(</>) and the entropy AH ((f)). 
In addition, we compute the angles between all pairs of vectors in the list and plot them as 
a histogram approximating the distribution function g(<p). 

B. Results for quantum kicked top 

For our numerical study of the kicked top, we use a twist parameter k = 3 and either spin 
J = 511.5 or spin J = 63.5. The calculations are done in the corresponding 512-dimensional 
or 64-dimensional even-parity subspace; the even-parity subspace is the eigenspace of the 
parity ( |3.4j ) with eigenvalue +1. Throughout this section, we consider only two initial states. 
The first one, the coherent state \9, <p) with 9 and (p given by Eq. (|3.6j ), is centered in a regular 
region of the classical dynamics; we refer to it as the regular initial state. The second one, 
referred to as the chaotic initial state, is the coherent state \9,(p) with 9 and (p given by 
Eq. ( |3.7| ); the chaotic initial state is centered in a chaotic region of the classical dynamics. 

Figure [I] shows results for spin J = 511.5 and a total number of 2 12 = 4 096 vectors after 
n — 12 perturbed steps. The interaction strength is g = 0.003, corresponding to perturbing 
rotations by an angle ±0.003 radians about the z axis (the perturbation strength in units of 
the coherent-state radius 2J" 1 / 2 is gJ 1 ^ 2 /2 = 0.034). For Fig. |l|(a), the chaotic initial state 
was used. The distribution of Hilbert-space angles, g(<p), is concentrated at large angles; 
i.e., most pairs of vectors are far apart from each other. The information AI needed to 
track a perturbed vector at resolution level is 12 bits at small angles, where each group 
contains only one vector. At <fi ~ 7r/16 the information suddenly drops to 11 bits, which 
is the information needed to specify one pair of vectors out of 2 11 pairs, the two vectors in 
each pair being generated by perturbation sequences that differ only at the first step. The 
sudden drop of the information to 10 bits at </> ~ 7r/8 similarly indicates the existence of 2 10 
quartets of vectors, generated by perturbation sequences differing only in the first two steps. 
Figure |I|(a) suggests that, apart from the organization into pairs and quartets, there is not 
much structure in the distribution of vectors for a chaotic initial state. The 2 10 quartets 
seem to be rather uniformly distributed in a = \2 AH2 ~\ = 46-dimensional Hilbert space, 
a supposition supported by a comparison with Fig. |^, where results for 2 13 = 8 192 random 
vectors in 45-dimensional Hilbert space are shown. 

The inset in Fig. |l](a) shows the approximate functional dependence of the information 
needed about the perturbation, A/ min , on the tolerable entropy AH to \, based on the data 
points AI(<f)) and AH {(f)). There is an initial sharp drop of the information, reflecting the 
grouping of the vectors into pairs and quartets. Then there is a roughly linear decrease of 
the information over a wide range of AH tQ \ values, followed by a final drop with increasing 
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slope down to zero at the maximum value of the tolerable entropy, AH to \ = AH$. The large 
slope of the curve near AH to \ = AH$ can be regarded as a signature of hypersensitivity to 
perturbation. The linear regime at intermediate values of AH to i is due to the finite size of 
the sample of vectors: in this regime the entropy AH r of the rth group is limited by log 2 N r , 
the logarithm of the number of vectors in the group. As is discussed in more detail in an 
upcoming publication |p8 |, finite-sample-size effects dominate the behavior of the curve as 



long as the number of vectors in a typical group is smaller than the effective dimension of 
Hilbert space. One therefore expects reduced finite-size effects for larger samples and smaller 
Hilbert-space dimensions. 

To investigate this prediction, we plot in Fig. [5] results for J = 63.5 and 2 14 = 16 384 
vectors after 14 perturbed steps, using an interaction strength g = 0.03 (gJ 1 ^ 2 /2 = 0.12). 
As expected, the inset shows that the large-slope behavior of the AJ min vs. AH to i curve 
extends to higher values of A/ min than in the previous case. Otherwise, Fig. is very similar 
to Fig. |I|(a). Apart from an obvious organization of the vectors into pairs, there is no 
apparent structure in the distribution of vectors. The 2 13 pairs of vectors are distributed 
quasi-randomly in a n& = 45 dimensional Hilbert space, as is confirmed by comparing Fig. |3| 
to the case of 2 13 random vectors in 45 dimensions shown in Fig. [2|. 

We limit the discussion of the regular initial state to spin J = 511.5 since, as pointed 
out in Sec. [Ill B|, only for a large dimension of Hilbert space is the regular initial state well 



localized inside the regular region around the elliptic fixed point (|3.6j ). Figure |l](b) shows 



data for 2 12 vectors after 12 perturbed steps in the regular case. The distribution of per- 
turbed vectors starting from the regular initial state is completely different from the chaotic 
initial condition of Fig. 0(a). The angle distribution g(<fi) is conspicuously non-random: it is 
concentrated at angles smaller than roughly 7r/4, and there is a regular structure of peaks 
and valleys. Accordingly, the information drops rapidly with the angle <fi. The number of 
explored dimensions is n& = 2, which agrees with results of Peres [§] that show that the quan- 
tum evolution in a regular region of the kicked top is essentially confined to a 2-dimensional 
subspace. The A/ m i n vs. AH to \ curve in the inset bears little resemblance to the chaotic 
case. Summarizing, one can say that, in the regular case, the vectors do not get far apart 
in Hilbert space, explore only few dimensions, and do not explore them randomly. 

To obtain better numerical evidence for hypersensitivity in the chaotic case and for the 
absence of it in the regular case would require much larger samples of vectors, a possibility 
that is ruled out by restrictions on computer memory and time. The hypothesis most 
strongly supported by our data is the random character of the distribution of vectors in 
the chaotic case. To close this section, we give a crude argument to show that randomness 
in the distribution of perturbed vectors implies hypersensitivity to perturbation. Consider 
N vectors distributed randomly in a .D-dimensional Hilbert space, where N ^> D ^> 1, 
so that AHg ~ log 2 -D. In this case it was shown in that AI((f>) ~ — D log 2 (sin 2 0) and 
AH ~ sin 2 01og 2 D ~ sin^Aif^ if AI(cf)) < log 2 N. The information is thus proportional 
to the dimension of Hilbert space, D, whereas the tolerable entropy is proportional to log 2 D. 
Eliminating <fi gives 

(4.9) 
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i.e., hypersensitivity to perturbation. 



V. QUANTUM MAPS ON THE UNIT SQUARE 
A. Quantized unit square 

In this section, we investigate perturbed quantum maps on the unit square. Our quanti- 



zation procedure closely follows p9f and [p0|| . To represent the unit square in D-dimensional 
Hilbert space, we start with unitary "displacement" operators U and V, which produce dis- 
placements in the "momentum" and "position" directions, respectively, and which obey the 
commutation relation 

UV = VUe , (5.1) 

where e D = 1. We choose e = e 2m / D . We further assume that V D = U D = r\ = ±1. The 
case r] = 1 corresponds to periodic boundary conditions (used in Sec. |VE|) , whereas the case 



7] = — 1 corresponds to anti-periodic boundary conditions (used in Sections |V C] and |V D| ). 
It follows |29|j30[| that the operators U and V can be written as 

U = e 2m « and V = e~ 2 ^ . (5.2) 

The "position" and "momentum" operators q andp both have eigenvalues (/3 + j)/D, where 
the integer j runs over a sequence of D integers and 

/*={?' ^ V = \ (5-3) 

i.e., T] = e 2m/3 . 

This quantization procedure fixes the eigenvalues of q and p to lie within some unit 
interval, but leaves us free to choose the interval. We use the freedom to write the eigenvalues 
of q and p as 

qj = P±l= p ., j = Do,...,D + D-l, (5.4) 

where the integer Do can be chosen freely. We denote by \qj) and \pj) the eigenvectors of q 
and p with these eigenvalues. In Sections |V C| and [VD| , we choose D = 0, corresponding to 
qj,Pj G [0, 1), whereas in Sec. |V b\ we choose D = —D/2, corresponding to qj,Pj G [— |, |). 



In the following, we represent Hilbert-space operators by matrices in the position basis 
For consistency of units, let the quantum scale on phase space be 2nTi — 1/D. A 
transformation between the position basis {|<?j)} and the momentum basis {|Pj)} is effected 
by the operator Gn, defined by the matrix elements 

(G D ) kj = (p k \q 3 ) = e ~ ip ^ /h = J— e -*«W+W+i)/D _ (5>5) 

v D 



Following [30|, we define a family of D 2 states \qj,Pk) {j, k = D , . . . , D + D — 1), which 
we call, in analogy to the continuous case, coherent states. The fiducial coherent state \qo,Po) 
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(assuming — D < D < 0) is denned as the eigenstate with the smallest eigenvalue of the 
operator [30] 

W = 2 - ^ = 2 - cos(27rg) - cos(27rp) . (5.6) 

The other coherent states are defined by applying the displacement operators to the fiducial 
state, i.e., 

\qj,Pk) = e~ injk/D U k V j \q ,p ) = e^ D V j U k \q ,p ) . (5.7) 



The coherent states can also be characterized as minimum- uncertainty states [ 30] . 
The coherent states are an overcomplete set, satisfying 

jjJ^llj'P^iQi'P^ = 1 = ( unit operator) . (5.8) 

Therefore they can be used to define a discrete Q function, representing an arbitrary state 
[■0) via 

Q\i>)(qj,Pk) = -^\{(ij,Pk\ip)\ 2 ■ (5.9) 

In terms of the Q function, the coherent state \qj,Pk) is represented by roughly circular 
contours centered at the point (qj,Pk), and it occupies an area of about 1/D pOfl . 

For our discussion of the different quantum maps below, we always assume the initial 
state \ipo) to be a coherent state. The state \ip n ) after n unperturbed steps is given by 

\ijn}=f n \^} } (5.10) 



where T = Tb for the baker's map (Sec. [V Q ), T = Tl for the lazy baker's map (Sec. [VD| 
and T = Ts for the sawtooth and cat maps (Sec. |V~E]). 



B. Perturbation operator 
In Sec. prO, we presented an example of an interaction with an environment that can 



be described as a stochastic perturbation of the unitary time evolution. Here we assume 
from the start that at each step a random perturbation operator is applied following the 
unperturbed evolution operator T. The perturbation operator is constructed by analogy 
to our previous work [|J on the classical baker's map, where the perturbation consisted 
in shifting the phase-space pattern by a small amount in the p direction. The quantum 
analogue of this perturbation would be an operator that shifts the Q function of a state by 
a small amount in the p direction. Due to the discrete nature of our quantum version of the 
unit square, however, a Q function can be shifted rigidly only if it is shifted by a discrete 
multiple of 1/D. 

In order to apply a perturbation that does not shift the momentum by a discrete multiple 
of 1/D, one might consider a perturbation operator P = e - 2m P a e n a l h = e -2mpa e 2mqDa 
ramps up the phase uniformly in the position basis, i.e., P^j = <5£;je _27n/3o e 27n<?jDa = 5/ CJ e 27r4: ' a . 
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When a is not a multiple of l/D, however, this perturbation can couple widely separated 
momentum eigenstates because of the abrupt phase change between j = D + D — 1 and 
j = D . To avoid this abrupt phase change, we use instead a perturbation operator II 
defined in the position basis by 

fr _U fcj e 2 ^ a , if Do < j < D + D/2, . . 

kj ~ \ 5 kj e 2 ^ 2D °+ D -^ a , if D + D/2 < 3 <D Q + D, ^ AAJ 

where a > is the magnitude of the momentum shift (Da is the magnitude of the momen- 
tum shift in units of the separation between momentum eigenstates) and where we restrict 
attention henceforth to Hilbert spaces with an even dimension D. Roughly speaking, the 
operator II shifts the momentum up by a in the left half of the unit square and down by a 
in the right half of the unit square. 

A perturbed time step consists of first applying the unperturbed time-evolution operator 
T, followed either by the perturbation operator II or by its inverse IT -1 , chosen randomly. 
After n perturbed time steps, the number of different perturbation sequences — or histories — 
is 2 n , all occurring with equal probability. 



C. Baker's map 

The classical baker's transformation maps the unit square < q,p < 1 onto itself 
according to 



1 (2 9 -l,i(p + l)), ifi<«<l. 

This corresponds to compressing the unit square in the p direction and stretching it in the q 
direction, while preserving the area, then cutting it vertically, and finally stacking the right 
part on top of the left part — in analogy to the way a baker kneads dough. 

There is no unique way to quantize a classical map. Here we adopt the quantized baker's 
map introduced by Balazs and Voros |31] and put in more symmetrical form by Saraceno 
30fl . We assume anti-periodic boundary conditions (77 = —1, f3 = |) and D = 0, which 



leads, as discussed above, to discrete position and momentum eigenvalues at half-integer 
values qj = (| + j)/D — Pj, j = 0, . . . , D — 1. The quantum baker's map is now defined by 
the matrix 

where D is even and the matrix elements are to be understood relative to the position basis. 

We simulate the perturbed quantum baker's map for two different Hilbert-space dimen- 
sions, D = 16 and D = 256, but using the same perturbation parameter a = 0.005 and the 
same initial coherent state |^o) = \1/2D,1/2D) (i.e., j = k = 0) in both cases. In units 
of the separation between p eigenstates, this perturbation is 16 times as big for D = 256 
(Da = 1.28) as for D — 16 (Da = 0.080). 

Figure |](a) shows the results for D = 16. Because the perturbation is so small, the 
vectors after 14 perturbed steps are not spread very far apart in Hilbert space, a typical 
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angle between vectors being = tt/8. The information decreases rapidly with resolution 
angle 0, and the entropy AH$ is less than 1 bit. Nevertheless, the data have several fea- 
tures characteristic of chaos: (i) no structure is apparent in the distribution function g(<p), 
indicating a random distribution of vectors within the small region explored by the vec- 
tors; (ii) even though the entropy AHs is very small, the number of explored dimensions is 
rid = 6; (iii) although the A/ min vs. AH to \ curve in the inset is qualitatively different from 
the chaotic cases in Figs. [I] and |3|, it shows nonetheless that in order to keep the entropy 
increase AH to \ below 10% of its maximum value AH$ requires at least 12 bits of information 
about the perturbation. These three features should be contrasted with the corresponding 
behavior for the regular initial state of the quantum kicked top [Fig. 0(b)]. 

The effects of the perturbation become more pronounced if the Hilbert-space dimension 
is increased to D = 256, as can be seen in Fig. ^(b), where results for 2 13 vectors after 13 
perturbed steps are plotted. Here, besides an organization into pairs of vectors a distance 
of about 37r/16 apart, the distribution of the vectors is similar to the distribution of ran- 
dom vectors in Fig. |2]. The vectors explore rid — HI Hilbert-space dimensions, and the 
AI m i n vs. AH to \ curve in the inset has a steep slope near AH to \ = AHs, a signature of 
hypersensitivity to perturbation. 



D. Lazy baker's map 



The lazy baker's transformation maps the unit square < q,p < 1 onto itself 
according to 

if < q < I, 

{q,p) i— { (l-p,q) , rt\<q<\, (5.14) 




if ~ < q < I- 



Here the unit square is cut into three vertical stripes, the left and right stripes are compressed 
in the p direction and stretched in the q direction by a factor of 3, the middle stripe is rotated 
by 90°, and finally the three parts are stacked on top of each other — somehow imitating the 
action of a lazy baker who neglects to knead the middle portion of the dough. 

To define a quantized lazy baker's map, the unit square is quantized in the same way as 
for the ordinary baker's map (rj = —1 and D = 0). The dimension D of Hilbert space is 



chosen to be a multiple of 6. The quantum-mechanical time-evolution operator is [32 










D/3 








Gd/3 



(5.15) 



where I d/3 is the identity matrix in D/3 dimensions. 

The quantum lazy baker's map was introduced [Q because it allows investigation of 
regular as well as chaotic evolution, while being almost as simple as the ordinary baker's 
map. As we did for the kicked top, we compare regular and chaotic evolution by choosing 
different initial vectors. We use a Hilbert-space dimension D = 216 — even multiples of 



powers of 3 reflect best the symmetries of the map j32] — and a perturbation parameter 
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a = 0.005 (Da = 1.08) to facilitate comparison with the results of the last section for the 
ordinary baker's map. 

For the chaotic case shown in Fig. |5J(a), we choose the initial state to be the coherent 
state \if) ) — 1 52. 5/216, 52.5/216) (i.e., j = k = 52). The Q function of this state, of diameter 
~ 1/ a/216 ~ 1/15, covers a region containing a large number of hyperbolic points separating 
regular regions in a fractal hierarchy ||32|| ; thus one expects chaotic behavior for this initial 
condition [52]. Indeed, Fig. |5J(a), similarly to Fig. f|(a), shows that after 13 steps the 2 13 



vectors are distributed quasi-randomly, aside from an organization into pairs of vectors. 
The vectors explore n d = 136 dimensions, and the A/ min vs. AH to \ curve has the shape that 
seems to be characteristic for chaos. 

The results for the regular case, shown in Fig. ||(b), are strikingly different. For this 
simulation, we choose the initial state to be the coherent state that is specified by j = k = 108 
(i.e., |V>o) = |108.5/216, 108.5/216)). This state, also of diameter 1/^216 ~ 1/15, is almost 
entirely confined to the central square of side length 1/3, which is simply rotated by 90° 
at each step. One expects regular behavior for this trivial time evolution, and Fig. |5](b) 
confirms this expectation. The peaks and valleys of the distribution function g((j)) indicate 
a regular spacing between the vectors, the vectors explore only = 2 dimensions, and the 
A-fmin vs. AH tol curve has the shape that we encountered in Fig. |I|(b), which seems to be 
characteristic for regular dynamics. 



E. Sawtooth and cat maps 



The classical sawtooth map |33| on the torus q,p & [— |, \) is given by 



(q, p) 1 — ► \ {q + p + Kq) mod l,(p + Kq) mod 1 



(5.16) 



where the sawtooth parameter K is an arbitrary real number and mod 1 is taken relative 
to the relevant interval [— 5, 5). For K = 0, the map is trivial and completely regular, 
describing a free rotor. For K > 0, the map is completely hyperbolic (chaotic). For K — 1, 
it becomes the Arnol'd cat map [34|. 

For the quantized version of the sawtooth maps [ 55155 1, we assume periodic boundary 
conditions (77 = 1, (3 = 0) and D = —D/2, which leads to discrete position and momentum 
eigenvalues at qj = pj = j/D, j = —D/2, . . . ,D/2 — 1, where D is even. The quantum 



evolution operator is given by the matrix elements p5| , p6 



SJjk 



\Ts\qk) 



— 27T /4 

£ e i-wKk 2 1 D e iir{j-k) 2 / D 



(5.17) 



For the simulation of the perturbed sawtooth map, we use the Hilbert-space dimension 
D = 256, a very small perturbation parameter a = 0.001 (Da = 0.256), and the initial 
coherent state l^o) — [ 0, 0) . We compare the distribution of 2 13 vectors obtained after 13 
perturbed steps for three different values of the sawtooth parameter. 

For Fig. H](a) we make a generic choice of the sawtooth parameter, K = 0.5, corresponding 
to chaos in the classical map. As in the behavior found for the 16-dimensional baker's map in 
Fig. [|a, the very small perturbation strength keeps the vectors from moving very far apart. 
Nevertheless, the vectors explore = 38 dimensions, and to keep the entropy increase 
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to 10% of its maximum value requires nearly 12 bits. Altogether, as for the baker's map 
in 16 dimensions, we come to the conclusion that the sawtooth map with K = 0.5 shows 
characteristics of hypersensitivity to perturbation. 

In contrast to the generic sawtooth map, consider the free rotor resulting from a sawtooth 
parameter K — 0, a case illustrated in Fig. |6](b). The free rotor shows entirely regular 
behavior. Although the perturbation manages to move some vectors far apart from each 
other, as is evident from the shape of the distribution function g(<f)), the vectors explore only 
n d = 2 dimensions of Hilbert space. The A/ min vs. AH to \ curve has the shape already found 
in the regular cases of the kicked top and the lazy baker's map. 

Finally, Fig. [7] shows results for the Arnol'd cat map, generated by a sawtooth parameter 
K — 1. Integer values of the sawtooth parameter K are non-generic in the sense that 
they give the map additional symmetries These additional symmetries might 

be responsible for the structure in the angle distribution g(<p), apparent in Fig. [7], which is 
absent in the generic case shown in Fig. ^(a). Another difference between the cat map and 
the generic sawtooth map, which might also be due to the additional symmetries, is the 
grouping into pairs, quartets, and possibly octets, revealed by the sudden drops in the AI 
vs. cf) curve in Fig. [7|, but which is absent in Fig. |6|(a). As already discussed for the kicked 
top in Sec. [IV BL this grouping shows up in the A/ min vs. AH to \ plot as an an initial drop of 
the information. Apart from this initial drop, however, the AJ min vs. AH to \ curve is similar 
to the generic case, indicating that the quantum cat map, like the generic sawtooth map, is 
hypersensitive to perturbation. 



VI. CONCLUSION 

Our aim in this paper is to explore, via numerical simulations, the notion of hypersen- 
sitivity to perturbation in quantum systems, trying to determine whether hypersensitivity 
might provide a criterion, as it does for classical systems, for distinguishing regular from 
chaotic quantum dynamics. What has been learned? Our simulations indicate that hyper- 
sensitivity is a valuable tool for investigating chaotic quantum dynamics. Our key finding 
is that if a chaotic quantum system is perturbed stochastically, the different perturbation 
histories tend to produce state vectors that are distributed randomly over a large number 
of Hilbert-space dimensions, whereas if a regular quantum system is subjected to a similar 
perturbation, the resulting state vectors explore only a few Hilbert-space dimensions and do 
not explore even those dimensions randomly. 

The nearly random distribution of vectors translates into a characteristic A/ min vs. AH to \ 
curve, evident in the insets of Figs. [j](a), |2], |3], [|(b), and §(a), which seems to be a property 
of quantum chaos. This characteristic curve is distinctly different from the A/ m i n vs. AH to \ 
curve produced by regular quantum dynamics [Figs. 0(b), 0(b), and ||(b)]. After an initial 
sharp drop due to perturbation histories that differ only at the first few steps, the A/ min 
vs. AH to \ curve for chaotic dynamics has a region of roughly linear decrease that extends over 
a broad range of values of AH to \, this linear decrease is a consequence of the finite sample of 
vectors produced by the perturbation (always a binary perturbation in our simulations) and 
thus cannot be said to be a property of the chaotic dynamics. The most important part of 
the A/ min vs. AH tol curve is the steep increase in A/ min as Aif to i is decreased just below its 
maximum value of AH$. This steep increase means the following: to decrease the entropy a 
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small amount below the entropy that comes from averaging over the stochastic perturbation 
requires a large amount of information about the perturbing environment, much larger than 
the entropy decrease. This rapid increase in information is what we call hypersensitivity to 
perturbation. It seems to be a feature of chaotic quantum dynamics. 

These considerations set the stage for further work: to understand in detail the steep 
rise in the A/ m i n vs. AH to \ curve and how it is cut off by a finite sample of vectors; to show 
that quantum hypersensitivity to perturbation is indeed a property of chaotic quantum 
dynamics and not an artifact of particular perturbations; to study the time dependence 
of quantum hypersensitivity and perhaps to extract from the time dependence a quantity 
like the Kolmogorov-Sinai entropy of classical chaotic dynamics; and, finally, to develop a 
general theory of quantum hypersensitivity to perturbation, which would allow one to relate 
it to other properties of chaotic quantum dynamics. 
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APPENDIX A: 

Equation ( |2.5| ) is true for any measurement on the environment that is described by a 
set of completely positive operations jF r ; an operation is a linear, trace-decreasing mapping 
that takes positive operators to positive operators and, hence, takes a density operator to 



an (unnormalized) post-measurement density operator JT7j. The joint state of the system 



and the environment after a measurement that yields results r is 

•7>(Ptotal) _ J2i AriPtotalAli 
Pr Pr 

The complete positivity of T r means that it can be written as the sum involving the operators 
A ri ||17|| ; since we are considering a measurement on the environment, these operators are 
environment operators. The probability to obtain measurement result r becomes 

p r = tr(V r (/3 t otai)) = tv(p toisl E r ) , 

where 

R - V 4 1 " A ■ 

i 

is an element of a POVM. Completeness of the POVM guarantees that 

^ E r = lLg = (environment unit operator). 

r 

The system state after a meaurement that yields result r is obtained by tracing out the 
environment, 
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a _ tr £ (^ r (ptotal)) _ tr £ (p totalJ B r ) 

Pr — — • 

Pr Pr 

The crucial step here uses the cyclic property of the trace to move the operators A ri to 
the right side of the trace; this step is allowed within an environment trace because these 
operators are, by assumption, environment operators. The important content of Eq. ( |A1|) 
is that even though there are many operations T r that are consistent with the POVM E r , 
the system state after a measurement on the environment depends only on the POVM. 
Equation (|2.5|) follows from summing over r and using the completeness of the POVM: 

Y^VrPr = trg(ptotal) = f>S ■ 
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FIGURES 



FIG. 1. Results characterizing the distribution of Hilbert-space vectors for the perturbed kicked 
top with J = 511.5, k = 3, and g = 0.003 after n = 12 time steps. The two main diagrams show, as 
a function of the angle <j>, the distribution g((p) of Hilbert-space angles (unnormalized, in arbitrary 
units), the average information AI((f>) to specify a vector at the resolution given by (f> (in bits), 
and the average conditional entropy AH ((f)) (in bits). See the text for a precise definition of 
these quantities. In the insets, AI(<ft) is plotted versus AH ((f)) using the same data points as 
in the main diagrams. The labeling of the inset axes is motivated by the fact that AI((p) is a 
good approximation to A/ m i n for AH to \ = AH ((f)). (a) Chaotic case, i.e., initial coherent state 
IV'o) = |$> V 9 ) centered in the chaotic region with 6 and tp given by Eq. 
vectors generated after 12 perturbed steps, (b) Regular case, i.e., initial coherent state centered at 
the elliptic fixed point given by Eq. ( |3.6| ). All 2 12 vectors generated after 12 perturbed steps. 

FIG. 2. The same quantities as in Fig. | plotted using a list of 2 13 random vectors in 45-di- 
mensional Hilbert space. 

FIG. 3. As Fig. |(a) (chaotic case), but using J = 63.5, g = 0.03, and all 2 14 vectors generated 
after 14 perturbed steps. 

FIG. 4. As Fig. |l|, but for the quantum baker's map with initial coherent state 
\ifj ) = \1/2D,1/2D). (a) D = 16, a = 0.005, all 2 14 vectors after 14 steps, (b) D = 256, 
a = 0.005, all 2 13 vectors after 13 steps. 

FIG. 5. As Fig. [l], but for the quantum lazy baker's map. D = 216, a = 0.005, all 2 13 vectors 
after 13 steps, (a) Chaotic case, initial coherent state \ifjo) = |52. 5/216, 52.5/216). (b) Regular 
case, initial coherent state \if)o) = 1 108. 5/216, 108.5/216). 

FIG. 6. As Fig. ||, but for the quantum sawtooth maps. Initial coherent state | V'o) = 1 0, 0) , 
D = 256, a = 0.001, all 2 13 vectors after 13 steps, (a) Chaotic sawtooth map, K = 0.5. (b) 
Regular free rotor, K = 0. 

FIG. 7. As Fig. [l], but for the quantum cat map, i.e., the quantum sawtooth map with K = 1. 
Initial coherent state \if)o) = |0,0), D = 256, a = 0.001, all 2 13 vectors after 13 steps. 
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